Spin-orbit coupling effect in (Ga,Mn)As films: anisotropic exchange interactions and 

magnetocrystalline anisotropy 
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The magneto-crystalline anisotropy (MCA) of (Ga,Mn)As films has been studied on the basis of 
ab-initio electronic structure theory by performing magnetic torque calculations. An appreciable 
contribution to the in-plane uniaxial anisotropy can be attributed to an extended region adjacent 
to the surface. Calculations of the exchange tensor allow to ascribe a significant part to the MCA 
to the exchange anisotropy, caused either by a tetragonal distortion of the lattice or by the presence 
of the surface or interface. 

PACS numbers: 75.50.Pp, 75.30.Gw, 73., 75.70.-i 



Diluted magnetic semiconductors (DMS) are a class 
of materials having attractiveproperties for spintronic 
applications (e.g. see review Many investigations 

in this field are focussed on the (Ga,Mn)As DMS sys- 
tem with 1 to 10% of Mn atoms which have promising 
features from a physical as well as technological point 
of view. The crucial role of valence states with respect 
to various magnetic properties of (Ga,Mn)As was dis- 
cussed in the literature by many authors [l|, Q. First 
of all, the valence band holes are responsible for ferro- 
magnetic (FM) order in the system mediating the ex- 
change interaction between well localized Mn magnetic 
moments. Spin-orbit coupling of the states at the top 
of valence band, being close to the Fermi level, leads 
to a rather strong cubic magnetocrystalline anisotropy 
(MCA) in bulk (Ga,Mn)As and to an in-plane biaxial 
MCA in the (Ga,Mn)As film on top of a GaAs sub- 
strate In the latter case the spin-orbit coupling 
(SOC) makes the valence states close to Ep sensitive to 
lattice distortions and in that way responsible for the in- 
planc MCA due to compressive strains originating from 
the lattice mismatch between the (Ga,Mn)As film and 
GaAs substrate [i-14|. As soon as the spin polariza- 
tion of the valence bands is rather small, the MCA in 
(Ga,Mn)As is discussed in terms of anisotropic exchange 
interactions of the Mn atoms j^-Q ■ The strength of the 
MCA depends on the hole concentration introduced by 
the Mn impurity atoms H, B 15, EBl as well as on the 
variation of equilibrium lattice parameter of (Ga,Mn)As 
increasing with the increase of Mn content and resulting 
in a larger lattice mismatch with the GaAs substrate. 

Numerous experimental results evidenced a transition 
from the bi-axial to the uni-axial in-plane anisotropy 
15|. So far, however, there is no consensus in the lit- 
erature concerning the origin of the in-plane uniaxial 
anisotropy. Although in some recent theoretical works 
the origin of the uniaxial in-plane anisotropy is attributed 
to a trigonal distortion caused by a uniaxial or shear 
strain within the film plane IH, [13 ; ttiis type of dis- 
tortion was not observed experimentally. Stacking fault 



defects in the (111) and (111) planes have been found re- 
cently in experiment [18| which could be responsible for 
breaking the equivalence of the [110] and [110] directions 
in the (Ga,Mn)As films. However, there is so far no ex- 
perimental evidence nor theoretical description showing 
that these stacking faults are responsible for the in-plane 
uniaxial anisotropy. 

In order to obtain a more detailed understanding of 
the subtle electronic effects which determine the MCA 
properties of (Ga,Mn)As films, we performed ab-initio 
electronic structure calculations for tctragonally dis- 
torted (Ga,Mn)As bulk as well as (Ga,Mn)As films de- 
posited on a GaAs substrate. The ab-initio calculations 
have been performed within the framework of the lo- 
cal spin density approximation (LSDA) of density func- 
tional theory (DFT) using the fully relativistic Korringa- 
Kohn-Rostoker (KKR) multiple scattering band struc- 
ture method [l^ 0] ■ For the treatment of the chemi- 
cal disorder in (Ga,Mn)As alloys we applied the coher- 
ent potential approximation (CPA). Moreover, for the 
bulk and surface calculations we used a regular /c-mesh 
of 63 X 63 X 63 points in the full 3D Brillouin Zone (BZ) 
and 63 x 63 points in the full 2D BZ, respectively. For the 
angular momentum expansion of the Green's function a 
cutoff of ^niax = 3 was applied. 

The study of the magneto-crystalline anisotropy 
(MCA) was performed by calculating the magnetic 
torque T^'^'^ = —dE{{ek})/dei x acting on the mag- 
netic moment rfii of the atomic site i, with a unit vector 
Bi — mi/\rhi\ pointing along the direction of the magne- 
tization M . The component of the magnetic torque with 
respect to the axis u 



n{e,<j>)^-dE{M{e,<j>))/de 
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was calculated from first-principles as described in j2l| . 
Here, the u vector specified by the angles 9 and </> (see 
Fig. [T^) lies within the surface plane and is perpendicu- 
lar to the direction of the magnetic moment cm- For an 
uniaxial anisotropy a special geometry can be used which 
gives a simple relationship between the magnetic torque 
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and the energy difference between the in-plane and out- 
of-plane magnetization directions. Setting 9 = 7r/4, the 
torque component gives the dependent energy dif- 
ference Tu{9 = 7r/4, (/)) = — E±. In the case of an 
in-plane anisotropy these values can also be used to eval- 
uate the anisotropy energy within the plane, comparing 
in particular the directions [110] and [110]. 

The exchange coupling tensor J_^^ used below for the 
discussions of the magnetic anisotropy in terms of the 
anisotropic Mn-Mn exchange interactions [H-Q was cal- 
culated as described in Ref. [l^]. Here, the effective co- 
efficients of the uniaxial MCA are represented by the fol- 
lowing expression [H, 01 : 



[001] 
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with Ki being the on-site MCA coefficients [23|. 

In order to study the strain-induced effect in the MCA 
of deposited (Ga,Mn) As films, we consider at first a bulk 
system with tetragonal distortion (avoiding surface and 
interface contributions) which is then characterized by 
the c/a ratio. Magnetic torque calculations simulat- 
ing the strain-induced effects in the alloy with 5% Mn 
yield a linear variation of the magnetic anisotropy energy, 
i?[ioo] — £-[001] J from -1-3.38 to —3.37 /ieV per unit cell for 
c/a ratio varying from 0.99 to 1.01, i.e. the magnetic easy 
axis changes from an out-of-planc to an in-planc orien- 
tation, which is in line with corresponding experimental 
data [25[. As the [100] and [010] directions are equiva- 
lent, this leads to the bi-axial in-plane MCA with [100] 
and [010] being easy magnetization directions and an in- 
plane anisotropy energy i?[ioo] ^ -£-[110] ~ — 0.1 yueV per 
unit cell. 

For a more detailed analysis of the relationship be- 
tween the MCA and anisotropy of Mn-Mn exchange in- 
teractions, calculations of the exchange coupling tensor 
elements J^"^ have been performed for (Ga,Mn)As with 
5% Mn both without any distortion as well as with a 
tetragonal distortion of c/a = 1.01. For an undistorted 
(Ga,Mn)As system we find that the sum EjX-^g" " ^ij^) 
(a, /3 = x,y, z) over all lattice sites in the expression given 
in Eq. ^ vanishes (see Fig. [T|d). This is a consequence 
of the system's symmetry, in spite of the fact that the in- 



dividual terms (J"" — Jf!^) with a ^ B are non-zero 
the presence of a tetragonal distortion along the z-axis, 
the symmetry properties within the xy plane (i.e. (001) 
plane) do not change. Therefore, summation over all lat- 
tice sites up to Rij = 5a (with lattice parameter a) shown 
in Fig. [TJ: gives ^jiJif' — JfJ') = 0. For more details. 

Fig. [TJl shows the differences J|^-^ — JfJ' for Rij taken 
along [100] and [010] directions (dashed lines). These 
values are finite and equal in magnitude, but they have 
an opposite sign and therefore cancel each other upon 
summation over all sites. However, due to the tetragonal 
distortion along z, {J^j^ — JfJ') for Rij taken along [001] 
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FIG. 1. a) Geometry for the torque calculations; b) J^^ — J^^ 
and J^^-jyy for bulk (Ga,Mn)As with 5% Mn, with tetrago- 
nal distortion c/a = 1.01; c) T^jiJ^f-J^ and Y^jiJ^-J^D 
over all lattice sites up to Rij < 5a; d) J"^"" — J^^ along [100] 
and [010] directions. 



and [010] directions are not equivalent (see Fig. [TJl, solid 
lines) and thus the sum J2ji'^if~'^if) O'^^r all lattice sites 
does not vanish anymore. The summation over all lattice 
sites up to Rij < 5a is shown in Fig. [TJ; which gives the 
contribution to the uniaxial MCA that originates from 
the exchange anisotropy being « 2.5 ficV. Because of the 
slow convergence of the sum with increasing distance, 
this gives only an approximation to the true contribu- 
tion due to the exchange anisotropy. Nevertheless, the 
value obtained in this way has the same order of magni- 
tude as the MAE obtained from our torque calculations 
leading to the conclusion that the exchange anisotropy 
has indeed a significant impact on the total MAE. 

Our present investigations of the in-plane uniaxial 
anisotropy have been performed for a 8 monolayer (ML) 
thick (Ga,Mn)As film deposited on a semi-infinite (001)- 
oriented GaAs substrate. In order to distinguish the 
anisotropy behaviour in the vicinity to the interface with 
GaAs as well as in the area adjacent to the surface we 
performed calculations for an uncovered (Ga,Mn)As film 
as well as one with two additional capping layers of An. 
Due to small amount of free charge carriers in (Ga,Mn)As 
the surface potential decays slowly into bulk leading to 
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a potential and a charge density gradient within an ex- 
tended region adjacent to the surface. The existence of 
such a potential gradient results in the breaking of the 
4-fold symmetry of the bulk (Ga,Mn)As system, making 
the [110] and [110] directions inequivalent (for the geom- 
etry used here this corresponds to the x and y directions, 
respectively) and leading effectively to a C2v symmetry 
not only within the few surface/interface layers but also 
in a rather extended subsurface regime. 

We discuss now the surface induced MCA in the 
film. Here we focus mainly on the MAE properties of 
a (Ga,Mn)As film with a clean Ga terminated surface 
deposited on GaAs(OOl). The results for the energy dif- 
ferences between different magnetization directions are 
^^[110] - -£^[001] = -80.56 /xeV and E^^^q^ - £^[001] = 
—32.96 /xeV per film unit cell (8 ML). This gives an uni- 
axial in-plane anisotropy with the energy difference of 



110] 



-47.6 /ieV per film unit cell. 



£^[110] - 

Fig. [2^ presents the layer resolved contributions to 
the i?[iio] — £'[001] a-iid iJ^ixg] — £'[001] values, indicated 
by open and filled symbols, respectively. The difference 
between these values characterizes the MCA within the 
plane. One should emphasize here that the contribu- 
tion to the MCA from the region close to the surface de- 
cays slowly into the bulk. Therefore the surface-induced 
anisotropy effect in the uniaxial in-plane MCA is deter- 
mined by a rather extended region adjacent to the sur- 
face and not just by two or three subsurface layers as 
it is often observed in metallic systems (e.g. 



261). The 



corresponding contribution to the energy of the uniaxial 
in-plane anisotropy exceeds by far the energy of the bi- 
axial in-plane anisotropy when normalized to the same 
volume (£'[100] — £'[110] ~ —0.1 /xeV per unit cell of the 
bulk system). Using these results the MCA of exper- 
imental (Ga,Mn)As films consisting of 7i + 8 monolay- 
ers can be modelled by combining the contribution of n 
bulk-like layers with the contribution of 8 layers of sur- 
face region. This gives two competing contributions to 
the MCA: a bi-axial in-plane anisotropy from bulk-like 
layers of (Ga,Mn)As with a tetragonal distortion and a 
uniaxial in-plane anisotropy from the area adjacent to 
the surface. Applying our obtained MAE values to a 
unit volume, one can get the MCA of the whole film 
including the surface region. Within our consideration, 
the coefficient of the in-plane bi-axial anisotropy does 
not depend on the film thickness L, while the the coef- 
ficient of the in-plane uni-axial anisotropy K2 recalcu- 
lated per unit volume should decrease with film thick- 
ness as 1/L. Thus, according to our numerical results, a 
rather strong uniaxial anisotropy should be observed in 
the case of very thin films, while the increase of the film 
thickness should lead to a competition of bi-axial and uni- 
axial anisotropics beginning with a certain film thickness. 
The contribution from the 'surface' region to the out-of- 
planc uniaxial anisotropy £^[iyo] ~ ^[001] decreases as well 



with the film thickness as 1/L. This results in a leading 
role of the in-plane anisotropy contribution caused by 
the tetragonal lattice distortion discussed above. Note 
that an increase of the Mn concentration results in an 
increase of the charge carriers in the film which again 
results in better screening of the surface potential. This 
can be seen in Fig. [2^, where the values £'[110] — £'[001] 
and £^[iYo] ~ £-[001] ^re shown as a function of the dis- 
tance from the surface for a (Ga,Mn)As film with 11% 
Mn. This increase in Mn concentration results in an in- 
plane MAE £'[110] — £'[iTo] ~ —20.8 ^eV per film unit 
cell, i.e. one obtains a smaller anisotropy energy when 
compared to the case of 5%Mn. 

Since the uniaxial MCA has its origin in an extended 
subsurface region one can expect that it is an intrin- 
sic property of the systems and should be observed not 
only in the case of a clean surface but also in the pres- 
ence of overlayers on the top of the (Ga,Mn)As film. 
Corresponding investigations have been performed for a 
(Ga,Mn)As film with 3 capping layers of Au on top of 
the (Ga,Mn)As film. The resulting layer resolved con- 
tribution to the MCA is shown in Fig. [^J). In spite 
of the differences in the MAE between the Au capped 
(Ga,Mn)As film and the case of uncovered film, the gen- 
eral trend in both cases is the same, i.e. one can clearly 
see that the difference in layer contributions to the MCE 
for different directions of magnetization, along [110] and 
[110], decays slowly with the distance from the surface or 
Au/(Ga,Mn)As interface, respectively. 

To investigate also the effect caused by a concentra- 
tion gradient along the surface normal within an un- 
covered (Ga,Mn)As film we dealt with a corresponding 
film where the Mn concentration varies from 5% at the 
(Ga,Mn)As/GaAs interface to 6.6% in the surface layer. 
As can be seen in Fig. [2^ such a gradient does not result 
in a noteworthy change in the MCA. 

To analyze in more detail the origin of the surface- 
induced in-plane uniaxial anisotropy the contribution of 
the exchange interaction anisotropy in the (Ga,Mn)As 
film was determined. Fig. [3] shows the difference J^ '^^ — 
Jfj calculated along the [110] and [110] directions within 
the film layers where the x and y axes are chosen along 
[110] and [110] directions, respectively. As discussed 
above, for bulk (Ga,Mn)As the variation of Jf^~Jf^ with 

distance \Rij \ is the same for Rij along the [110] and [110] 
however with different sign. This behaviour is more or 
less the same for Mn atoms next to the (Ga,Mn)As/GaAs 
interface (see Fig. [3^). For Mn in the fourth layer (chos- 
ing the surface layer as the first layer), however, the situ- 
ation is changed indicating a pronounced modification of 
the anisotropic exchange coupling due to the broken sym- 
metry. As a result, the sum over lattice sites in Eq. ([2]) 
does not vanish which leads to a contribution to the in- 
plane uniaxial anisotropy. Fig. [Sb shows the correspond- 
ing results obtained by summing the terms (Jf,-^ — JiF) 
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FIG. 2. a) Layer resolved contributions to the MCA energy in 
the uncovered 8ML (Ga,Mn)As film with 5 at.% Mn (circles) 
and ll(triangles) at.% Mn, as well as with Mn content varying 
from 5 at.% at the interface to 6.6 at.% in surface layer {gradx, 
squares), for two directions of magnetization: M||[110] and 
M||[110]; b) layer resolved contributions to the MCA energy 
in the 8ML (Ga,Mn)As film with 5 at.% Mn, with 3 capping 
layers of Au. 

over all lattice sites j within a sphere of radius 2.9 a with 
i taken within the layers 1 — 8 in the (Ga,Mn) As film. As 
one can see, the anisotropy of the exchange interaction 
gives indeed a substantial contribution to the anisotropy 
energy i5[iio] — £'[iTo] ^'^^ ^^'^ layers 3 to 8. For the first 
two film layers, i.e. surface and subsurface layers the two 
curves strongly deviate reflecting the dominating on-site 
contribution to the MCA [23l.[24|. 

In summary our results show that the tetragonal dis- 
tortion (caused by a compressive strain due to lattice 
mismatch of (Ga,Mn)As and GaAs lattices) is respon- 
sible for the bi-axial in-plane anisotropy that is in line 
with the interpretation given in previous investigations. 
A strong uniaxial in-plane MCA was found in (Ga,Mn)As 
film in the area adjacent to the surface or an interface. 
We conclude that this is a result of the slow decay of the 
surface potential gradient due to the small amount of free 
charge carriers. The contribution to the uniaxial in-plane 
anisotropy decays rather slowly into the bulk and is not 
restricted to only a few surface layers. Moreover, a signif- 
icant contribution responsible for the MCA in the films is 
caused by the anisotropic Mn-Mn exchange interactions 
mediated by holes in the valence band of (Ga,Mn)As. 
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3. Variation of {JiJ^ 
of Mn atoms taken in [110] and [110] directions in the 
(Ga,Mn)As/GaAs film: a) bulk vs (Ga,Mn)As/GaAs inter- 
face and b) bulk vs (Surf.- 3)-layer ; c) Layer-resolved sum 
— "^jiJif — Jfj^) calculated within the sphere of radius 2.9a 
in comparison with the MCA energy ij'^^^l — ijl^ioi evaluated 
by magnetic torque calculations for the 8ML (Ga,Mn)As film 
with 5 at.% Mn. 
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